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Abstract. Medical ultrasound scanners are typically calibrated to the soft tissue 
average of 1540 m s“^. In regions of different sound speed, for example organs and 
tumours, the B-mode image then becomes a distortion of the true tissue cross-section, 
due to the misrepresentation of length and refraction. We develop a ray model to 
predict the geometrical distortion of a general elliptical object with atypical speed of 
sound, and quantify the ensuing image distortion as a function of the key parameters, 
including the speed of sound mismatch, the object size and its elongation. Our findings 
show that the distortion of area can be significant, even for relatively small speed of 
sound mismatches. Our findings are verified by ultrasound imaging of a test object. 
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1. Introduction 

Medical ultrasound imaging is routinely used to evaluate length, area and volume as 
part of clinical diagnostic assessment, for example in sizing of organs and tumours, 
and assessing the age, weight and normality of fetuses [T]. The images are obtained 
through ultrasound time-of-flight measurements, converted into the spatial domain via 
an assumed speed of sound, typically taken to be the soft tissue average of 1540 m s“h 
However, bodily soft tissue is acoustically inhomogeneous and its speed of sound can 
vary by several percent across different regions [ 2 ]. For example, typical speeds of sound 
in fat, muscle and liver are 1440, 1580 and 1590 m s“^, respectively [3]. The image 
is then geometrically distorted from the actual tissue profile, where the variable speed 
of sound is not accounted for. This, in turn, leads to systematic inaccuracies in the 
evaluation of length, area and volume of anatomical structures. 

To illustrate this effect, consider an ambient medium with speed of sound Vi and 
an ultrasound scanner which assumes the same value. An object with atypical speed of 
sound V 2 7 ^ Vi is embedded within the medium and a B-mode image taken of its cross- 
section. Then, the distance within the object will be misrepresented in the image and the 
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refraction at the object boundary will not be accounted for, both of which will contribute 
to the geometric distortion of the object’s image. While this effect has been predicted 
previously HEllellTllllE!, a quantitative understanding is lacking. Such understanding 
is important in establishing errors in clinical measurements, and devising corrective 
strategies. Take, for example, the evaluation of the volume of organs and tumours now 
possible through three-dimensional ultrasound imaging [IDIIIIIII2]. This capability has 
diagnostic importance, particularly in charting the evolution of a disease, and has been 
demonstrated on a variety of structures including the bladder da m nans], prostate 
HJI, kidney ra. thyroid ra and tumours m- The 3D image is typically reconstructed 
from 2D images [ 211122 ] ; then, if the object in question has atypical speed of sound, 
the above 2D distortion will introduce errors in the calculated area of each scanning 
plane, and thus the total calculated volume. To date, this distortion remains largely 
unexplored as a source of error in area and volume measurements. 

In this work we make a quantitative study of the geometric distortion arising in 
2D ultrasound images of a general elliptical object with atypical speed of sound. We 
develop a ray model to predict this distortion (presented in Section , applicable across 
the four main multielement transducers currently in medical use (linear, curved, phased 
and vector arrays). Our aim is not to predict the full intensity distribution of the 
ultrasound images; to do so would require a sophisticated model of the acoustic field 
to take into account wave effects (such as diffraction and interference) and the nature 
of the beam line (such as its finite width and focussing). Our aim instead is to predict 
and understand the underlying mapping of the object boundary to the image boundary; 
as such, the ray approach is well-suited to this aim. Its simplicity also lends to future 
development of strategies to correct the distortion. Comparison to ultrasound images of 
phantoms (described in Section]^ confirm the success of the ray model in predicting the 
shape of the image. We analyse and quantify the distortion as a function of the speed 
of sound ratio V 2 /V 1 , the object’s shape and size, and the array type. Our results show 
that a small speed of sound mismatch can lead to significant distortions: for example, 
a 5% mismatch in the speeds of sound can lead to around 15% error in the area of a 
circular object, and this is sensitive to the size and shape of the object. 

2. Theory 

2.1. Ray Model 

We consider an elliptical (two-dimensional) object in the x-z plane, with uniform speed 
of sound ^ 2 , embedded in an ambient medium with uniform speed of sound v^. The 
object has semi-widths b and c. The set-up is depicted in figure 3(a). The scanner is 
assumed to be correctly calibrated to the ambient speed of sound vi. We will develop 
the model for the case of a curved transducer array; later, we will generalize to the other 
main array types. In such an array, beam lines can be considered to spread out from 
a virtual source within the (curved) transducer, emanating normal to the transducer 
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Zob(x) 


Figure 1. (a) Schematic of the curved transducer array and elliptical object (speed 
of sound V 2 ) embedded within an ambient medium (speed of sound vi). For V 2 ^ ui, 
refraction and misrepresentation of length within the object cause a point C on the 
far boundary to appear at C. (b) A close-up of the real ray (black) and image ray 
(grey), and the angles involved, at the point B. 


surface j23]. The origin is taken to be the virtual source. For generality we allow for a 


lateral offset a between the transducer and the ellipse. The front face of the transducer 


lies a distance d from the source. The object boundary is then represented by the 
elliptical line Zoh(x) which satisfies the relation 



( 1 ) 


We model the distortion of the object by treating the beam lines as rays. This 
approximation is valid in the far-field of the source and for objects much larger than the 
acoustic wavelength. Also, where the beam is incident on a boundary, the ray approach 
ignores mode conversion to interface waves that might occur at certain crucial angles. 
Nonetheless, through comparison to experimental images, we will see that the ray model 
accounts for the geometric distortion in the image. 

A ray is emitted from the transducer surface at point A with coordinates {xa,za) 
and is incident on the front face of the object at point B with coordinates ( 0 : 3 , zb)- This 
follows the straight line path zab{x){x) = xcoIOa, where 6a is the angle of the emitted 
ray relative to the vertical 


( 2 ) 


9 a = arctan —. 

Zb 


At B, the angle between the normal of the boundary and the vertical, On, can calculated 
from the gradient of the boundary at B via 



c^{xb - a) 


( 3 ) 


dx (xb^zb) 


lA{d + do + c - Zb)' 


The angle of incidence at B (relative to the normal) is then di = + ^a- 
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At the front face, a component of the wave will be back-scattered to the transducer, 
registering the presence of this point in the image. Since this ray travels only through 
the ambient medium, for which the scanner is assumed to be calibrated, the image of 
this point corresponds to its real location. Of more interest to us is the component of 
the wave which is refracted through the boundary (the components of the wave that 
undergo specular reflection or scattering away from the transducer are not registered by 
the transducer and so we ignore them). This ray is incident on the object’s back face 
at point (7, with coordinates {xc,zc), and represented by the straight line z^cix), 

ZBcix) = (x- xb) cot 9b + Zb. (4) 


6b is the angle from the vertical direction, given by 9b = 9i — 9n, where 9^ is the angle of 
refraction. To obtain an expression for 9^-, we rearrange Snell’s law, sindi/ sin^r = Vi/v 2 , 
to give 


9t: = arcsin 


— sin(6»A + 6'n) 
A 


( 5 ) 


To establish {xc^zc), we seek the intersections of the refracted ray zbc{x) with the 
boundary 2 ;ob(^)- Rearranging equation @ for X and inserting into equation 0 gives 
the relation 

[{z - Zb) tan 9b + xb - of [z - c - d - dof _ . . 

P ^ 

Multiplying this out leads to a quadratic equation in and thus two solutions for ; 2 . The 
larger of these corresponds to point C (and the smaller to point B). Hence we obtain 
the 2 ;-coordinate of C to be 


= Zb 


( 7 ) 


+2i 


c^{xb — a) tan 9b — b^d' — zbc^ tan^ 9b 
c2 tan^ 9 b + 


n 2 


c 262 — b‘^d'^ — (? [xb 


a 


Zb tan 6 *b] 


c2 tan^ 9b + b‘^ 


where d' = d + do + c. xq then follows from equation ([^ . 

Having established the path to (7, we proceed to locate its image, C'. The paths 

AB and BC have lengths Lab == a/x| -|- 2:3 and Lbc = \J{xc — xb)^ + {zq — zbY, 
respectively. Assuming that the ray returns to the transducer along the same path, the 
total travel time is 


_ 2Lab 2Lbc 

L “1“ . 


( 8 ) 

Vi V2 

The transducer, calibrated to the ambient speed of sound Vi, perceives the ray as the 
unrefracted path, i.e. along zp^b{x)^ with the same travel time t. Hence the path AC 
has length 

Vl 


Lac — — Lab H— -Lbc- 

2 V2 

The position of C then follows as 

(xc', Zc) = {Lac sin^A, Lac cos 9a). 


( 9 ) 


( 10 ) 
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For vi 7 ^ V 2 the image point differs from the real point, partly due to refraction of the 
ray at the near face of the object, and secondly due to the rescaling of length due to 
the speed of sound misrepresentation. For V 2 < ui, the image is extended from the 
transducer (due to the longer travel time within the object), while for V 2 > vi the image 
is moved closer to the transducer (due to the shorter travel time within the object). 


2.2. Generalization to Other Array Types 

There are 3 other main types of multielement transducer currently in use in medical 
imaging: the linear, phased and vector arrays. The ray structure for these array types 
can be found elsewhere m- Below we summarize how the ray model can be extended 
to these types. 


Vector array: In this array, beam lines spread out from a virtual source located within 
the (flat) transducer. Since we assume that vi = uq, no virtual refraction takes place 
at the transducer surface, and so the only modification is in the position at which 
the ray emerges from the surface, accounted for in the model by the coordinates 
{xa,Za). 

Phased array: In a phased array beam lines spread out from a virtual source at the 
surface of the transducer, accounted for by setting do = 0. 


Linear array: In a linear array, beam lines emerge from the transducer surface, i.e. 
do = 0, normal to the transducer face, i.e. 0^ = 0. Note that the image position 
(10) then reduces to 


{xc',Zc') = (0, Z/AcO) 


( 11 ) 


i.e. the image point appears directly above the source. 


2.3. Area Evaluation 


The image boundary will, in general, be irregularly-shaped. To evaluate the enclosed 
area we exploit Green’s theorem [23j. Green’s theorem states that, for a region R bound 
by a closed path F, which is continuous and with no crossing points, then, 




da: dz = 


{u dx -\- V dz} , 


r 


( 12 ) 


da: dz, which 


where u{x,z) and v{x,z) are arbitrary functions, and the path of integration along F 
is taken to be counter-clockwise. Taking the choice of functions u{x,z) = —zl2 and 

v{x,z) = x/2 then the left-hand side of equation (12) reduces to JJ 

corresponds to the area of the region R. It then follows that we can express the area A 
of this region via a line integral around its boundary F via 


r 


{—z da; -|- a; d.^). 


(13) 
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Figure 2. (a) Illustration of the experimental set-up. A curved array transducer, 
clamped in place, images (from above) the central cross-section of a spherical bauble. 
The bauble is filled with liquid of speed of sound U 2 , and the ambient speed of sound 
is vi. (b) Corresponding 5-mode image for V 2 = vi = 1540 m s“^, both matched to 
the scanner. The expected object boundary is shown by the dashed green line. 


We numerically implement this line integral around the image boundary (discretised as 
a series of points) to obtain the image area. We quantify the image distortion through 
the area distortion A-JA,., where Ai and A,, are the areas of the image and object. 

3. Materials and Methods 

A hollow spherial aluminium bauble, filled with “liquid 2” and sealed with a latex 
condom, was suspended (by a wire frame) in a bath of “liquid 1”. Its central circular 
cross-section was imaged from above by an ultrasound scanner (Philips HDI 5000) 
and curved array (Philips C5-2 40R, with do — 40 mm). The set-up is depicted in 
figure [^a) . The bauble and condom are sufficiently thin to have negligible impact on 
our analysis. The transducer is aligned manually with the equator of the object. The 
overall object radius was measured with calipers as (3.99 ± 0.01) cm and the transducer- 
object distance d determined using the scanner’s electronic calipers. The experiments 
were run at ambient temperature and the precise bath temperature recorded by digital 
thermometer. 

For the liquids we employed pure water, pure methanol and an ethanol-water 
mixture (with 7.25% w/w ethanol). These have room temperature speeds of sound 
of 1482 m s“^ [25], 1100 m s“^ [26] and 1540 m s“^ [27], respectively (the latter being 
matched to the tissue average and scanner-assumed value). The bath temperature varied 
in the range 18-23°C and we corrected the speeds of sound used in the model using the 
established temperature-speed relations [251 ESI EZ]. 

As an initial validation we filled both the object and bath with ethanol-water 
mixture, such that the system has uniform speed of sound throughout, V 2 = vi = 1540 
m s“^. The 5-mode, shown in figure [^a), shows a circular image boundary, confirming 
the absence of geometric distortion, and matches well to the expected circular shape of 
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the object (green dashed line). Note the intensity distortion around the boundary: the 
poles of the object are brightest as a result of specular reflection. 

Three variations of this system were considered. 

V 2 < vii Taking ethanol-water mixture as liquid 1 and pure water as liquid 2 embodies 
a system with V 2 < Vi (v 2 ~ 0.96ni). 

V 2 Vii Taking the ethanol-water mixture as liquid 1 and now methanol as liquid 2 
embodies a system with V 2 -C vi {y 2 ~ 0.71ni). 

V 2 > Vi'. Taking pure water as liquid 1 and the ethanol-water mixture as liquid 2 creates 
a system with V 2 > vi {v 2 ~ l.Odui). 

Note that, for the first two cases, the speed of sound of the ambient liquid is 
matched to that assumed by the scanner, as per the assumptions of our ray model. This 
is not true for the third case, and so the ray model must be extended to account for the 
distortion of the near-face of the object; this is outlined in [Appendix A , 

4. Results: Circular Objects 

We begin by using the ray model to map out the area distortion of circular objects, first 
for a simplest case of a linear array and then for our primary case of a curved array. 
The latter will be validated against experimental images. Unless otherwise stated, we 
assume the nominal parameters r = 4 and d = 2 . Note that length is subsequently 
presented in arbitrary units. 

4 . 1 . Linear Array 

In figure a) we show the ray model prediction for V 2 /v\ = 0.95; the typical ray 
behaviour is illustrated by three sets of rays. The rays propagate vertically downwards 
from the transducer surface (the a:-axis, z = 0 ), impinging the near-face (the upper half) 
of the object boundary. Our assumption that the scanner is matched to the ambient 
speed of sound means that the near face is undistorted, that is, the image boundary 
(solid green line) and the object boundary (dotted black line) overlap here. 

At the near-face the real rays (red lines) are refracted inwards; the exception is the 
central ray which continues vertically due to its zero angle of incidence, Oi. These rays 
travel slower within the object since V 2 < Vi. The image rays (dashed blue lines) deviate 
in two ways: they are unrefracted at the near-face, and they appear elongated within 
the object. The latter effect is best resolved along x = 0: here the image ray is simply 
extended by a distance 2 r{yi/v 2 — 1) ~ 0.42. Off-axis, the rays are refracted inwards 
(towards the normal), with the refraction being small across most of the object since 6 ^ is 
small (with the deviation being hard to see by eye). Towards the equator of the object, 
however, d\ grows rapidly, causing a significant increase in refraction and the real-ray 
path length in the object. This has two notable effects as a: —?• ±r: firstly, the increasing 
inwards refraction leads to the appearance of a gap in the image boundary, and secondly. 
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Figure 3. (a) Image (solid green line) of a circular object (dotted black line) with 
V 2 = 0.95ui based on a linear array. Three example sets of real rays (dashed blue lines) 
and their corresponding image rays (solid red lines) are overlaid, (b) Same as above 
but for V 2 = l.OSui. (c) Area distortion A \/as a function of the speed of sound ratio 
V 2 IV 1 , for the linear array (green dashed line) and curved array (blue solid line). 


the divergence of the internal path length leads to the formation of prominent “flicks” 
in the image boundary [5] . 

Now consider the case V 2 /V 1 = 1.05, shown in flgure|^b). Since V 2 > Vi the real 
rays travel faster than the image rays within the object, and so the far image boundary 
becomes shifted towards the transducer. Also, the real rays become refracted outwards 
(away from the normal). This removes the appearance of the above flicks for a: ~ ±r. 
The image boundary again features a gap, although its origin is different to that above: 
here it is caused when the rays impinge the near face at a sufficiently large angle of 
incidence that they undergo total internal reflection, and so do not reach the far face. 

For the above cases the area distortion is A-JA^. = 1.093 and 0.907, respectively. In 
other words, a 5% speed of sound difference has led to ~ 10% error in the area. We 
have also explored the wider dependencies of the area distortion on the key parameters. 
The area distortion scales approximately linearly with V 2 lv\, as shown in flgure|^c). 
As d, the transducer-object distance, is varied, the image area A, is unaffected, since 
all rays approach the object in the same direction. Moreover, the image area scales in 
proportion to the size of the circular object, such that the area distortion ratio Ai/Aj. is 
independent of r, the radius of the object. 

4 . 2 . Curved Array 

The corresponding images for a curved array are shown in figure |^a-b). We see 
similar qualitative effects as for the linear array: the shifting of the image boundary 
in front/behind the object boundary, gaps in the image boundary, and the appearance 
of flicks for V 2 < Vi. However, the ray structure and image shape are modified further 
due to the angular spread of the emitted rays. The near-face of the object has reduced 
extent, such that a greater proportion of the boundary is distorted. The incident rays 



















Geometric distortion of area in medical ultrasound images 


9 



Figure 4. (a)-(b) As figure]^ for a curved array (with d + do = 2). (c) Area distortion 
Ai/Ar of the circular object under curved array imaging as a function of the ratio r/d. 
Various values of V 2 /V 1 are considered ranging from 0.9 (top line) to 1.1 (lowest line) 
in steps of 0.05. The values V 2 /V 1 = 0.95 and 1.05 are highlighted by solid lines. 


now approach the tangent of the object fore of the equator, which causes a corresponding 
shift in the gaps and flicks in the boundary towards the transducer. For the V 2 < vi 
case we find A-JA,. — 1.147; for V 2 > vi we find A^jAy. = 0.859. In other words, a 5% 
difference in sound speeds for the curved array leads to over 14% deviation in area. 

As for the linear array, the variation of the area distortion with V 2 /V 1 is 
approximately linear [figure |^c)], but with greater deviation than for the linear array. 
The dependence on r and d is more complicated for the curved array since these 
parameters control the angular spread of the incident rays. Figure |^c) shows that 
the change in area distortion (relative to the undistorted value Ay/Ay = 1) increases 
with r/d; this is due to the increasing angular spread of rays as r/d increases, leading 
to greater refractive effects and increasing the proportion of the boundary which is 
distorted. 

4 . 3 . Experimental Validation 

Experimental B-mode images (with curved array) of the circular test object are shown 
in figurealongside the ray model predictions, for following 3 speed of sound scenarios. 

(a) V 2 < vi'. Here we consider vi = 1540 m s“^ (ethanol-water mixture) and V 2 = 1490 
m s“^ (water), giving V 2 /V 1 = 0.97. No geometric distortion is evident at the near 
face but the far face is shifted away from the transducer, in good agreement with 
the ray model (dashed green line). 

(b) V 2 -C ui: Now taking V 2 = 1100 m s“^ (methanol) gives the more extreme case of 
^ 2 /'^! = 0.71. The far-face of the image is highly distorted, in good agreement with 
the ray model. 

(c) V2 > vp. Lastly, we take Vi = 1490 m s“^ (water) and V2 = 1540 m s“^ (ethanol- 
water mixture), giving V 2 /V 1 = 1.03. The ambient liquid is no longer matched to 
the scanner calibration, causing the near face to become slightly shifted towards 
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Figure 5. B-mode images of the circular object of speed V 2 (imaged from above in 
a liquid with speed Ui), with the ray model superimposed (green dashed line), (a) 
V 2 = 1490 m s~^ and vi = 1540 m s“^, giving V 2 /V 1 = 0.968. (b) V 2 = 1100 m s“^ 
and vi = 1540 m s“^, giving V 2 /V 1 = 0.71. (c) vi = 1490 m s“^ and V 2 = 1540 m s“^, 
giving V 2 /V 1 = 1.034. 



Figure 6. Curved array images (solid green line) of an elongated elliptical object 
(dotted black line) with 5=1 and c = 4 (d + do = 2). Shown are the results for (a) 
V 2 = 0.95ui and (b) V 2 = 1.05ui. 


the far face. Again, the overall image shape agrees well with the ray model (in this 
case in based on the extended theory in Appendix A). 

These results validate the use of the ray model in predicting the geometric distortion of 
the object boundary. 

5. Results: Elliptical Objects 

We extend our theoretical analysis to the more general case of elliptical objects, which 
allows us to explore the effect of elongation. For brevity, we focus on the curved array. 
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We first consider an elliptical object that is elongated along ; 2 , with nominal 
parameters 6 = 1 and c = 4. The image distortion, presented in figure for both 
(a) V 2 = 0.95ni and (b) V 2 = l.OSui, shows similar qualitative structure for the circular 
object [figure]^ but with enhanced distortion. For V 2 < vi, the flicks are even more 
pronounced than for the circular case, caused by the fast roll-off of the object boundary 
with X (which causes the angle of incidence to increase more quickly with x). The 
area distortions are Ai/Aj. = 1.58 (for V 2 /V 1 < 1) and 0.68 (for V 2 /V 1 > 1); this is 
vastly greater than for circular objects in Section m demonstrating the sensitive role 
of elongation. 

In contrast, for a flattened ellipse [figure (a)-(b), nominal parameters 6 = 4 and 
c = 1] the image distortion is reduced, with negligible flicks appearing in the image 
due to the slower roll-off of the boundary with x. Also, the narrowness of the object 
in z means that the rays have little distance to accumulate any significant distortion. 
The area distortions here are Ai/Aj. = 1.06 (for V 2 /V 1 < 1) and 0.95 (for V 2 /V 1 > 1), 
confirming the reduced distortion for flattened objects. 

To further quantify the dependence on elongation, figure (c) plots the area 
distortion Ai/Aj- as a function of ellipticity h/c. We clearly see that the distortion 
diverges for objects which are increasingly elongated (along z), and converges towards 




Figure 7. (a)-(b) Same as figurej^but for a flattened elliptical object with 6 = 4 and 
c = 1. (c) Area distortion Ai/Ar versus ellipticity b/c (with 6 fixed to unity) under 
curved array imaging. Various values of V 2 /V 1 are considered ranging from 0.9 (top 
line) to 1.1 (lowest line) in steps of 0.05. 
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the undistorted result for increasingly flattened geometries. 

6. Discussion and Conclusions 

We have developed a ray model that describes the geometric distortion of two- 
dimensional ultrasound images (based on linear, curved, phased or vector arrays) of 
circular and elliptical objects due to speed of sound discrepancies. The geometric 
distortion arises from the refraction of the sound and the misrepresentation of length 
within the object. By comparing to ultrasound images of a test object we have validated 
that the ray model successfully captures the underlying geometric distortion of the 
image. 

The geometric distortion of ultrasound images of circular objects has been examined 
previously PE!, including by ray models PHj. Our analysis extends this previous body 
of work by quantifying the distortion of area, and thereby providing new insight into 
the errors in the ultrasonic evaluation of area and volume of anatomical features. Our 
work also extends previous analyses by considering elliptical objects, which allows us to 
explore the role of elongation. 

When V 2 < Vi, the image is larger than the object, and vice versa for V 2 > Vi. The 
ray model predicts the emergence of gaps in the image boundary, an effect which may 
contribute to the phenomenon of acoustic shadowing PPEllll. For V 2 < Vi the image 
also features distinctive flicks at the extreme angular positions, as noted previously [5|. 
The deviation of the image area from the true value is greater for the “angled” array 
systems (phased, curved and vector arrays) than the linear array, due to the enhanced 
role of refraction. The deviation also depends on the shape and size of the object, 
in general increasing with the lateral size of the object and its elongation along the 
ultrasound axis. 

One might expect that the area distortion follows the fractional mismatch in the 
speed of sound [7|. However, we typically find the distortion to be considerably larger. 
Take, for example, the circular object case considered earlier: there a 5% mismatch in 
the speeds of sound led to a 9% error in the area for a linear array and over 14% error 
for a curved array. The area distortion has contributions from refraction of the rays at 
the near surface and the miscalculation of the length of the ray within the object. We 
can artificially remove the effect of refraction from the ray model by setting O,- — d\. In 
doing so we find for a circular object that these effects roughly contribute to half of the 
area distortion (i.e., refraction accounts for 40-50% of the distortion); more generally 
this proportion will depend on the object shape. 

As well as geometric distortions, ultrasound images undergo well-known intensity 
distortions, arising from wave behaviour and the true beam line structure [9|. In 
future it would be insightful to consolidate these two distortion contributions by direct 
comparison between ray predictions and more sophisticated finite-element models of the 
full acoustic field. 

Our work has considered objects with known areas and derived their image. Of 
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more importance in the clinic is the inverse problem, that is, to obtain the true area 
given an image, and developments in this direction could promote algorithms to correct 
the distortion in clinical measurements. A further extension of this work would be 
to investigate the volumetric distortion of 3D objects with atypical speed of sound. 
Due to the increased dimensionality of this situation, one would expect the distortion 
to be further amplified compared to the 2D case. With ultrasound-based volume 
measurements now claiming to have errors of a few percent [isiiniin], these distortions 
may become the significant contributor to volumetric errors. 

Appendix A. Ray Model for Vi ^ Vq 

For simplicity our main ray model considered the ambient speed of sound vi to be 
matched to that of the scanner Uq, thereby negating the misrepresentation of distance 
as the ray travels through the ambient liquid. Here we extend the ray model to account 
for V\ 7 ^ Vq with a curved array. 


( 0 , 0 ) 



Figure Al. Schematic of the transducer-object scenario, featuring an elliptical object 
(speed of sound ^ 2 ) within an ambient medium (speed of sound ui). The transducer 
is calibrated to a speed of sound vq. 


We consider a transducer element on the surface of the curved array, point A 
Since the ray is emitted normal to the surface, no refraction effects occur 


[figure Al 


at the transducer boundary, despite vi ^ vq. The path of the real ray, which refracts 
into the object at point B on the near face and reaches the far boundary at point C, 
is provided by the same theory as before. The difference, however, now arises in the 
position of point B and its image point B’. We add a misregistration to all of the points 
by introducing the two relations for the misplacement of x and ; 2 , namely: 

Vq 


AA = I — _ 1 
vi 


+ z'^ — do 

sin6(A; AZ = ( — — 1 ) 

y/x"^ + z'^ — do 

- 

Vd / 

- 


cos 9 


A, 
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which are similar to those in Ref. j23]. Thus the position of the image point B’ is 


(xb',zb') — {xb + AX, Zb + AZ). 

The total length of the virtual ray, the equivalent of equation i§, 

LaC' = H- 

Vi V2 

and the image point C’ is located at 


IS 


(xcA^cO = {Lac' sin0A,^AC' cos^a). 


Thus the distortion of the near-face and the far-face is completely determined. 
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